Chirality in Bare and Passivated Gold Nanoclusters 
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Chiral structures have been found as the lowest-energy isomers of bare (Au28 and AU55) and thiol- 
passivated (Au2s(SCH3)i6 and Au3s(SCH3)24) gold nanoclusters. The degree of chirality existing 
in the chiral clusters was calculated using the Hausdorff chirality measure. We found that the index 
of chirality is higher in the passivated clusters and decreases with the cluster size. These results 
are consistent with the observed chiroptical activity recently reported for glutahione-passivated gold 
nanoclusters, and provide theoretical support for the existence of chirality in these novel compounds. 



PACS Numbers: 36.40. Cg,36.40.Mr,61.46.+w 

Detailed knowledge of the lattice structure, shape, 
morphology, surface structure, and bonding of bare and 
passivated gold clusters is fundamental to predict and un- 
derstand their electronic, optical, and other physical and 
chemical properties. This information is essential to op- 
timizing their utilization as novel nanocatalysts [GJ, and 
as building-blocks of new molecular nanostructured ma- 
terials 0, with potential applications in nanoelectronics 
H and biological diagnostics Q . An effective theoretical 
approach to determine gold cluster structures is to com- 
bine genetic algorithms and many-body potentials (to 
perform global structural optimizations), and first prin- 
ciples density functional theory (to confirm the energy 
ordering of the local minima) . Using this procedure we 
recently found many topologically interesting disordered 
gold nanoclusters with energy near or below the lowest- 
energy ordered isomer [^-^| . The structures of these clus- 
ters showed low spatial symmetry or no symmetry at 
all, opening the possibility of having distinct electronic 
and optical properties in such systems. In other stud- 
ies on passivated gold nanoclusters P,[l0|, we also found 
that the effect of a methylthiol monolayer (24 SCH3 
molecules) on a truncated-octahedron (with fee geome- 
try) AU38 cluster is strong enough to produce a dramatic 
distortion on the gold cluster, resulting in a disordered 
geometry for the most stable Au38(SCH 3 ) 2 4 passivated 
cluster. 

Although the calculated structure factors of the dis- 
ordered gold clusters were in qualitative agreement with 
the data obtained from X-ray powder diffraction on ex- 
perimental samples ||f?j], the direct confirmation of the 
existence of bare and thiol-passivated gold nanoclusters 
with low or no spatial symmetry had not been possi- 
ble due to the lack of enough experimental resolution for 
clusters in the size range of 1-2 nm JT^-|l|. Neverthe- 
less, in a recent study using circular dichroism, Shaaff 
and Whetten (SW) jl4| found a strong optical activ- 
ity in the metal-based electronic transitions (across the 
near-infrared, visible and near ultraviolet regions) of size- 
separated glutathione-passivated gold clusters in the size 



range of 20-40 Au atoms. SW pointed out that the most 
plausible interpretation of these results is that the struc- 
ture of the metal-cluster core of the gold-glutathione clus- 
ter compounds would be inherently chiral. Moreover, 
since the most abundant cluster in the experimental sam- 
ples corresponds to the passivated cluster Au2s(SG)i6, 
where SG denote the gluthatione adsorbate, SW pro- 
posed a chiral structure, with T point group, for the Au 2 s 
cluster 01. 



In earlier studies [|5| — |l 0|| , we have reported several struc- 
tural, vibrational and electronic properties of the disor- 
dered gold nanoclusters. However, although an initial 
quantification of the amount and type of local disorder 
present in amorphous-like structures of AU55 was ob- 



tained 1 15 1, to our knowledge, no attempt to theoretically 
investigate the existence of chirality in gold nanoclusters 
has been done before. In the light of the results ob- 
tained by SW regarding the intense chiroptical activity of 
size-separated glutathione-passivated gold nanoclusters, 
it is of interest to determine if the most-stable disordered 
structures we have found for bare and passivated gold 
nanoclusters are chiral, and to quantify their degree of 
chirality. In this letter, we present a structural analysis 
that shows the existence of chirality in these nanostruc- 
tures. This is done by implementing a method, based 
on the Hausdorff chirality measure p6| , Er| , to quantify 
the index of chirality of clusters. Our results show that 
the metal cores of the passivated clusters are more chiral 
than the bare clusters and that the degree of chirality de- 
creases with the cluster size. These results provide new 
insights on the effect of passivating agents on the bare 
metal cluster as a mechanism to generate chiral nanos- 
tructures with novel and perhaps unexpected properties. 

Our initial step was to determine the lowest-energy 
structures of the bare and passivated gold nanoclusters. 
To this end we first performed global optimizations using 
the Gupta n-body potential and a genetic-symbiotic al- 
gorithm for bare gold nanoclusters of different sizes (Ref. 
[0 provides details on this methodology) . Through this 
procedure we obtained the distribution in energy of the 
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most stable isomers for each cluster size. From this dis- 
tribution, we selected representative isomers like those 
with the lowest-energy, and those isomers that were con- 
sidered good candidates to be the lowest-energy minima 
based on the existence of well known symmetric struc- 
tures for certain cluster sizes like the truncated octahe- 
dron, icosahedron, decahedron, etc. The representative 
isomers were then locally reoptimized through uncon- 
strained relaxations using forces calculated from density 
functional theory (DFT) in the local-density (LDA) and 
generalized-gradient (GGA) approximations. The DFT 
calculations were done using the Siesta code |l8| with 
scalar- relativistic norm-conserving pseudopotentials [IE] 
and a double-C basis set of numerical atomic orbitals [IS] 
(see Refs. §JjJ for additional details on these calcula- 
tions). 

Here, we report results for the Au28 cluster since the 
lowest-energy structures for the Au„, n=38,55,75, clus- 
ters have been published elsewhere ■ Our results for 
the Au28 cluster show a similar pattern as was found for 
the larger sizes: the distribution of low-energy isomers 
shows a set of nearly degenerate in energy disordered 
structures, where the lowest-energy one corresponds to a 
highly distorted geometry. The T structure proposed by 
SW pi was obtained as a local minimum of the potential 
energy surface with higher energy. The top panel of Fig. 
1 shows both cluster geometries and the distribution of 
distances from the cluster center of mass. These struc- 
tures mainly differ in the presence of a 3-atom core in the 
lowest-energy disordered isomer whereas the T geometry 
contains a tetrahedral core. The DFT-LDA-GGA results 
confirm the higher stability of our disordered structure as 
compared with the T structure: we obtained an energy 
difference of -0.885 eV (LDA) and -0.638 eV (GGA) be- 
tween both isomers. A similar trend, showing the higher 
stability of disordered isomers with respect to ordered 
ones, has been obtained for larger cluster sizes: Au n , 
n=38,55,75 The physical origin of the higher sta- 

bility of amorphous-like structures in bare gold nanoclus- 
ters has been studied by analyzing the metallic bonding 
and its effect on the local stress of the clusters. It results 
that the amorphization is a mechanism for strain relief 
that lowers the cluster energy, being more prominent in 
the case of gold due to the short-range of the n-body 
metallic interaction 0S,Eq]. 

The lowest-energy structures of the thiol-passivated 
Au28(SCH 3 )i 6 and Au38(SCH 3 ) 2 4 clusters were obtained 
by performing local relaxations, using the forces calcu- 
lated from the DFT-LDA-GGA first principles method, 
and starting from different cluster-monolayer configura- 
tions. These include the lowest-energy ordered and dis- 
ordered bare Au28 and AU38 cluster geometries obtained 
by the procedure described above, with the methylthi- 
ols molecules placed on different adsorption sites (top, 
bridge, hollow), as well as on random positions over the 
metal cluster surface. Figure 2 (bottom panel) shows a 



highly distorted passivated cluster that corresponds to 
the lowest-energy structure of the Au 2 8(SCH 3 )i 6 cluster, 
obtained by the relaxation of the disordered bare Au28 
cluster with the thiols placed close to 3-atom hollow sites. 
From the distribution of radial distances to the cluster- 
center of mass shown in Fig. 2 (bottom panel), the degree 
of distortion and the ill-defined nature of the gold-thiol 
interface is evident. A similar strongly distorted geome- 
try was found for the lowest-energy structure of the larger 
passivated Au38(SCH 3 ) 2 4 cluster [p[[lO[]. The main driv- 
ing force producing the cluster distortion is the strong 
gold-sulphur interaction. An analysis of the metal-ligand 
interaction in the gold passivated clusters will be pub- 
lished elsewhere |pi]| . Meanwhile, recent calculations on 
smaller Au-thiolate compounds have confirmed that the 
strong covalent (directional) gold-sulphur interaction is 
mainly-responsible for the metal cluster distortion |^2| . 

The principal objective of this work is to determine 
if the disordered lowest-energy structures of the bare 
Au„, n— 28,38,55,75, and passivated Au28(SCH 3 )i6 and 
Au38(SCfi3)24 clusters are chiral, and to quantify the 
degree of chirality existing in them. This information is 
relevant for a proper interpretation of the circular dichro- 
ism measurements performed by SW on glutathione- 
passivated gold clusters |l4|] . Since chirality is a geomet- 
rical property of the system, independent of its chemical 
and physical manifestations, it is possible to quantify chi- 
rality without reference to experimental measurements, 
but using the inherent structural symmetry of the clus- 
ters. Although in recent years several approaches have 
been developed to measure chirality [ p3|p^ ], the Haus- 
dorff chirality measure, proposed by Buda and Mislow 
(BM) ]l6| ), has emerged as the general method of choice 
for the quantification of chirality jlTj]. Within this ap- 
proach, the degree of chirality is found by calculating the 
maximum overlap between the actual molecular struc- 
ture and its mirror image, using the Hausdorff distance 
p5~f between their sets of atomic coordinates. By rotating 
and translating one structure with respect to the other, 
the optimal overlap can be calculated. This Hausdorff 
chirality measure (HCM) is a continuous and similarity 
invariant function of the molecular shape, and is zero 
only if the molecule is achiral. The advantage of this 
approach is that its numerical implementation for large 
clusters containing n atoms in a three-dimensional space 
is straightforward, as was discussed by BM Our 
numerical procedure was tested by calculating the degree 
of chirality of the three tetrahedra shapes considered by 
BM, obtaining the same values as reported by these au- 
thors fl|. 

We calculate the HCM for the lowest-energy structures 
of the bare and thiol-passivated gold nanoclusters us- 
ing their relaxed Cartesian coordinates measured with 
respect to the cluster center of mass. Through an in- 
version operation, the coordinates of the mirror image 
clusters were obtained. The HCM was obtained by cal- 
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culating the maximum overlap between a given cluster 
and its mirror image. This corresponds to the mini- 
mum value of the Hausdorff distance between the sets 
of atomic coordinates of both structures. In order to ob- 
tain the minimum Hausdorff distance, the mirror cluster 
was translated and rotated around the original cluster 
in the three-dimensional space generating different con- 
figurations. For each of them, the Hausdorff distance 
with respect to the original cluster was calculated. The 
minimum of these values, normalized by the largest inter- 
atomic distance in the cluster, corresponds to the HCM. 
To check the reliablity of this methodology, the minimiza- 
tion of the Hausdorff distance was performed using the 
Broyden-Fletcher-Goldfarb-Shano method as suggested 
in Refs. (I^Jl^], a conjugate gradient optimization, and a 
direct evaluation of the Hausdorff distance taking varia- 
tions of 7r/360 for each rotational angle. The HCM values 
agree for the three methods up to the third significant 
digit. 

Our results show that the HCM index of chirality for 
AU38 and AU75 is zero. These values are expected since 
the lowest-energy structures corresponding to these sizes 
have one plane of symmetry and therefore are achiral. 
In fact, these results constitute an additional test of our 
numerical implementation for calculating HCM in large 
clusters. The HCM index of chirality for the lowest- 
energy structure and the T isomer proposed by SW for 
Au28 are displayed in Table I, together with the value ob- 
tained for the lowest-energy structure of AU55 (see Fig. 
I in Ref. JlB]]). These results show that the most stable 
structure of Au28, obtained in this work, is less chiral 
than the one proposed by SW, and that the index of 
chirality varies slightly with the cluster size from 28 to 
55 atoms. For the calculation of the HCM of the thiol- 
passivated gold nanoclusters, we only consider the coor- 
dinates of the Au atoms since we are interested in testing 
the suggestion of SW who proposed that the metal clus- 
ter core would be inherently chiral []l4| . The HCM values 
for the metal core of the passivated Au2s(SCH3)i6 and 
Au38(SCH3)24 clusters are also shown in Table I. These 
results show that both metal cluster cores are more chiral 
than the bare Au28 and AU38 clusters, and that the index 
of chirality for the passivated clusters decreases with the 
cluster size. It is noteworthy that the HCM values ob- 
tained for the bare and passivated chiral Au clusters are 
smaller than the values of the chiral tetrahedral shapes 
(HCM = 0.221, 0.252, 0.255 for D 2 , C 2 , and d tetrahe- 
dra, respectively) obtained by BM |f6) . In order to com- 
pare the HCM of Au clusters with a physical system, we 
calculated the index of chirality of the chiral D2-C76 and 
D2-C84, and the achiral Ih-Cea fullcrcncs, obtaining the 
values 0.109, 0.102, and 0.000, respectively. These results 
support the reliability of our methodology since there is 
experimental evidence on the existence of chirality in the 
-D2-C76 fullerene |^6j. This comparison also provides a 
useful reference to classify cluster and molecules in terms 



of their indexes of chirality, and indicates that the bare 
Au28 an d Au 55 clusters are as chiral as the .D2-C76, but 
the passivated Au28(SCH 3 )i 6 and Au38(SCH 3 ) 2 4 clusters 
are slightly more chiral than such fullerene. 

The above results on the existence of chirality in thiol- 
passivated Au nanoclusters are consistent with the chi- 
roptical activity measured on size-separated glutathionc- 
passivated gold nanoclusters in the size-range 20-40 Au 
atoms since it was shown that the optical properties are 
insensitive to the tail of the adsorbate thiol fll4j ]. Also, in 
the experimental samples it was found that the optical 
activity is weaker as the cluster size increases, and it is 
lost for larger clusters and for samples with mixed cluster 
sizes It remains to be verified if indeed there is a 

correlation between the value of the index of chirality and 
the strength of the intensity or other spectral features of 
the optical activity signal in the chiral samples. Another 
interesting prediction from the present calculations, to 
be confirmed experimentally, is related with the smaller 
index of chirality existing on bare clusters as compared 
to those for the metal-core of passivated clusters. Our 
results imply that the effect of the passivating monolayer 
is strong enough to distort a bare cluster geometry, pro- 
ducing chiral metal-cores that give rise to the intense 
chiroptical activity. This effect could change an achiral 
cluster into a chiral one, like in the AU38 case, or increase 
the index of chirality in an already chiral structure, like 
the Au28 cluster. 

The main contribution of this theoretical work has 
been to show the existence of, and quantify chirality in 
the lowest-energy structures of gold nanoclusters using 
the Hausdorff chirality measure. These results provide 
useful insights for a proper interpretation of the experi- 
mental results that report optical activity on passivated 
gold nanoclusters 0]. Our results also establish that 
passivating with thiol monolayers provides a mechanism 
to induce or increase chirality in bare Au clusters in the 
size-range of 20-40 atoms. This new information is rel- 
evant for the growth and sample preparation of differ- 
ent chiral metal nanoclusters. Theoretical precedents 
on metallic chiral nanostructures would include the NI28 
cluster |^], and the Al and Pb helical nanowires [B8[ , 
however the present results and those obtained by SW 
provide both theoretical and experimental evidence on 
the existence of chirality in gold nanoclusters. Electronic, 
optical, thermal, and other physical and chemical prop- 
erties of chiral gold nanoclusters are currently being cal- 
culated and will be reported in forthcoming publications. 
At present, it is expected that novel and interesting prop- 
erties emerge from the chiral character of metal clusters 
that could be useful for new nanotechnological applica- 
tions. 
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TABLE I. Index of chirality of the lowest-energy structures 
of bare (first column) and thiol passivated (third column) gold 
nanoclusters.The second column shows that all the ordered 
isomers are achiral except the T isomer of Au28- The index 
of chirality was calculated using the Hausdorff distance be- 
tween the sets of atomic coordinates of the cluster and its 
mirror image. The cluster geometry was obtained through 
structural optimizations using a many-body potential, genetic 
algorithms and density functional theory. 



Bare 




Passivated 


Au28 (disordered) 
0.106 


Au 28 (T) 
0.129 


Au28(SCH 3 )i 6 (disordered) 
0.160 


AU 38 (C a ) 

0.000 


Au 38 (O,0 
0.000 


Au38(SCHa)24 (disordered) 
0.121 


Auss (disordered) 
0.117 


Au 5 5 (h) 

0.000 




AU75 (Cs) 

0.000 


Au 75 (D 5h ) 
0.000 





FIGURE CAPTIONS 

FIG. 1. Top panel: Distances of the gold atoms from the 
center of mass and cluster geometry for the lowest- 
energy disordered (closed diamonds and inset (a)) 
and the T (open circles and inset (b)) isomers of the 
bare Au28 cluster. Bottom panel: Distances of the 
gold (open diamonds) and sulphur (stars) atoms 
from the center of mass for the lowest-energy struc- 
ture of the thiol-passivated Au28(SCH 3 )i 6 cluster. 
The closed diamonds denote the same gold atom 
distances as in the top panel. They are included to 
show at the same scale, the degree of distortion and 
expansion of the gold metal cluster core upon passi- 
vation. The insets show the geometries of the metal 
core and the passivated cluster. Sulphur atoms are 
depicted as darker spheres. 
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